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ABSTRACT 

We investigate the interactions of high-redshift galaxy outflows with low-mass virialized clouds of 
primordial composition. While atomic cooling allows star formation in objects with virial temperatures 
above 10'*K, "minihalos" with virial temperatures below this threshold are generally unable to form 
stars by themselves. However, the large population of high-redshift starburst galaxies may have induced 
widespread star formation in neighboring minihalos, via shocks that caused intense cooling both through 
nonequilibrium H2 formation and metal-line emission. Using a simple analytic model, we show that the 
resulting star clusters naturally reproduce three key features of the observed population of halo globular 
clusters (GCs). First, the lO^K maximum virial temperature directly corresponds to the ~ IO^Mq upper 
limit on the stellar mass of such clusters, a feature that can not be explained by any GC destruction 
mechanism. Secondly, the momentum imparted in such interactions is sufficient to strip the gas from 
its associated dark matter halo, explaining why GCs do not reside in the dark matter potential wells 
that are ubiquitous in galaxies. Finally, the mixing of ejected metals into the primordial gas provides 
a straightforward mechanism to explain the —0.1 dex homogeneity of stellar metallicities within a 
given GC, while at the same time allowing for a large spread in metallicity between different clusters. To 
study the possibility of such "fine grained" mixing in detail, we use a simple one-dimensional numerical 
model of turbulence transport to simulate mixing in cloud-outflow interactions. We find that as the 
shock shears across the side of the cloud, Kelvin-Helmholtz instabilities arise, which cause turbulent 
mixing of enriched material into ^ 20% of the cloud. Such estimates ignore the likely presence of large- 
scale vortices, however, which would further enhance turbulence generation. Thus the global nature of 
mixing in these interactions is multidimensional, and quantitative predictions must await more detailed 
numerical studies. 



1. INTRODUCTION 

Most of the Milky Way's ^ 150 globular star clusters are 
old. Their ages have been reliably determined from main 
sequence turnoff points, with the firm result that the ma- 
jority formed 10-13 Gyr ago and thus are coeval with the 
oldest stars in the Galaxy (Krauss & Chaboyer 2003, and 
references cited therein). Furthermore, the existence of 
similar systems in many other galaxies suggests that these 
objects represent an essential aspect of the epoch of galaxy 
formation. Seminal papers by Peebles & Dicke (1968), 
Searle & Zinn (1978), and Fall & Rees (1985) spurred nu- 
merous lines of investigation into this relationship and pro- 
duced scenarios in which globular clusters (hereafter GCs) 
formed either before, during, or shortly after the develop- 
ment of the galaxies that now host them. Unfortunately, 
improved theory and simulations, and especially, new ob- 
servations have made most of the earlier scenarios unten- 
able, or at least unlikely. Hence, the question of what 
triggers the formation of globular clusters has yet to be 
satisfactorily answered. 

One key development for this subject was the work of 
Zinn (1985), who established that two distinct populations 
of Galactic GCs exist: (1) relatively metal-rich clusters, 
with [Fe/H] > -0.8, which share spatial and kinematic 
properties with the Galaxy's thick disk or bulge; and (2) 
relatively metal-poor clusters, with -2.5 < [Fe/H] < -0.8, 



which appear to be part of the halo. Further work has in- 
dicated a similar dichotomy in other galaxies {e.g. Ashman 
& Bird 1993; Forbes, Brodie, & Huchra 1997; Beasley et 
al. 2000; Larsen et al. 2001), which suggests two different 
modes of GC formation. In the metal-rich case, forma- 
tion is likely to be an ongoing process. External galaxies, 
often in the midst of strong tidal interactions, have been 
observed to host OB associations with various properties 
expected of young GCs - they are compact, bright, bluish, 
and are estimated to have masses > lO^M© [e.g. Whit- 
more & Schweizer 1995; Schweizer et al. 1996). Further- 
more, it has recently been realized that galaxies can add 
to their system of globulars by stripping GCs from dwarf 
satellite galaxies, as seems to be occurring now in con- 
nection with the Milky Way's Saggitarius dwarf (Ibata et 
al. 2001; Yoon & Lee 2002). In low-metaUicity GCs, only 
ages > 10 Gyr are seen, and characteristics appear to be 
more uniform across galaxy types and sizes (see however, 
Strader, Brodie, & Forbes 2004). While the reader is re- 
ferred to the recent monograph by Ashman & Zepf (1998) 
for a more detailed comparison between these populations, 
our focus in this work will be to develop a formation model 
only for halo globular clusters. 

There are three essential features of these objects that 
any formation scenario must explain or accommodate. 
The first of these is the remarkable chemical homogene- 
ity of iron-peak elements exhibited by stars within a given 
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cluster: typically, the dispersion in [Fe/H] is less than 0.1 
dex {e.g. Suntzeff 1993). This is particularly extraordi- 
nary in the case of the oldest clusters, in which the star- 
forming gas must have attained its metallicity before the 
universe was a billion years old. This condition divides cur- 
rent formation theories for halo GCs into 'pre-enrichment' 
and 'self-enrichment' types. In the former scenario {e.g. 
Elmegreen & Efremov 1997; Bromm & Clarke 2002), GCs 
were formed out of gas that had already been homoge- 
neously enriched by a previous generation of supcrnovac 
(SNe). Here the key question is just exactly what that 
population was, why it played only a secondary role in 
the formation history of the GC, and how it could have 
enriched this material on very short time scales. In the 
self-enrichment picture, on the other hand, the protoclus- 
ter cloud was enriched by one or more supernova events 
occurring within it {e.g. Brown, Burkert, & Truran 1995; 
Cen 2001; Nakasato, Mori, & Nomoto 2000; Beasley et al. 
2003; Li & Burstein 2003). In this case, the key prob- 
lems arc that this sclf-cnrichmcnt must occur extremely 
homogeneously and that the kinetic energy corresponding 
to these SNe can be enough to unbind the gaseous proto- 
cluster (Peng & Weisheit 1991). Note however that the 
distribution of metallicities among different halo globular 
clusters is quite large, and approximates a Gaussian with a 
mean of [Fe/H] ~ —1.59 and a dispersion of [Fe/H] ~ 0.34 
(Ashman & Zepf 1998). 

A second constraint comes from observations of the tails 
of stars that are in the process of being stripped from GCs 
by the tidal field of the Galaxy (Irwin & Hatzidimitriou 
1993; Grillmair et al. 1995). If globular clusters contained 
substantial halos of dark matter, this increased gravita- 
tional potential would have a large impact on these tidal 
losses. In fact, no evidence of dark matter halo suppres- 
sion is seen, placing a strong upper limit of 2.5 on the ratio 
of total mass to stellar mass of these objects (Moore 1996; 
see however Maschenko & Sills 2004). This is in direct 
contrast with galaxies, which exhibit total mass to stellar 
mass ratios ~ 100 {e.g. Marinoni & Hudson 2002). Thus 
the mechanism that formed GCs is likely to be fundamen- 
tally different from the dark-matter driven collapse that is 
believed to govern the formation of galaxies, as it is much 
easier to separate dark matter from gas than from point- 
like bodies such as stars (see, however, Bromm & Clarke 
2002). 

The final constraint on GC formation is related to their 
mass distribution, which is well described as a Gaussian 
in logio(M*/M0) with a mean {M^) ^ lO^M© and a dis- 
persion of 0.5 {e.g. Armandroff 1989). This issue is most 
lucidly illustrated by the classic "survival triangle" in the 
mass-radius plane of globular clusters, which is bounded 
by the long-term disruption processes that act on Galactic 
GCs (Fall & Rees 1985). Here the minimum radius as a 
function of mass is determined by mechanical evaporation 
{e.g. Spitzer and Thuan 1972) while the maximum radius 
as a function of mass is bounded by shocking that occurs 
when a cluster passes through the high-density Galactic 
disk {e.g. Ostriker, Spitzer, & Chevalier 1972). While both 
these constraints seem to provide good agreement with the 
observed GC minimum stellar masses and sizes, dynami- 
cal friction, which is the only Galactic mechanism limit- 
ing the maximum stellar mass, only operates at masses 
> IO'^-^Mq. Thus it appears that the ~ lO^M© upper 



mass cut of globulars is not set by any known destruction 
mechanism {e.g. Gunn 1980; Caputo & Castcllani 1984), 
but rather represents an intrinsic property of the popula- 
tion of gaseous proto-clusters (Peng & Weisheit 1991). 

In the high-redshift universe, nature has provided us 
with just such a population. Because atomic line cooling 
is only effective at temperatures > lO^K, collapsed clouds 
of material with virial temperatures below this threshold 
must radiate energy through dust and molecular line emis- 
sion. While the levels of H2 left over from recombination 
arc sufficient to cool gas in the earliest structures {e.g. 
Abel, Bryan, Norman 2002; Bromm, Coppi, & Larson 
2002) the resulting 11.20-13.6 eV photon emission from 
the stars in these objects {e.g. Haiman, Rees, & Loeb 
1997; Ciardi et al. 2000) is likely to have quickly disso- 
ciated this primordial H2. Thus a generic prediction of 
current structure formation models is a large population 
of Tvir < lO-^K, Mgas < Kf'-^MQ virialized clouds of gas 
and dark matter that are unable to form stars until they 
interact with other objects. In fact, it was the similarity 
between these so-called "minihalos" and the globular clus- 
ter population that led Cen (2001) to propose that con- 
vergent ionization fronts during cosmological reionization 
might be able to transform these objects into GCs. Yet 
such a scenario begs the question of how globulars were 
enriched with metals, and it is not at all clear that the 
first galaxy-minihalo interactions were radiative. 

Observations of high-redshift starburst galaxies have un- 
covered large numbers of outflows, both in optical and in- 
frared measurements at 3 ^ 2; ^ 4 (Pettini et al. 2001) and 
in optical observations of lensed galaxies at 4 <J z ^ 6.5 
(Frye, Broadhurst, & Benitez 2002; Hu et al. 2002). In ad- 
dition, the classic picture of reionization is of a two stage 
process, which begins with individual sources ionizating 
their immediate surroundings and ends in a rapid "over- 
lap" phase in which neighboring H II region join together, 
quickly ionizing the remaining neutral regions {e.g. Gnedin 
2000). Thus it is possible that minihalos that were close to 
starbursting galaxies during the first stage of reionization 
could have been impacted by outflows prior to reioniza- 
tion fronts, depending on which was first able to escape 
from high-redshift starbursts. While numerical simula- 
tions have yet to yeild a definitive answer to this question, 
there are suggestions that the I-fronts in ^ lO^M© star- 
bursts may be D-type at small radii and slow to expand 
(Kitayama et al. 2004), and that shell material swept up by 
outflows is effective at trapping ionizing radiation (Fujita 
et al. 2003). And while eventually shell fragmentation al- 
lowed such ionizating radiation to escape, even then, large 
regions may have remained "shadowed" by the fragments. 

It is possible, therefore, that a large number of minihalos 
were first impacted by outflows. Unlike radiation fronts, 
which typically boil away the minihalo gas {e.g. Shapiro et 
al. 2004), shock interactions will result in intense cooling 
through both nonequilibrium formation of H2 (Mac Low & 
ShuU 1986; Shapiro & Kang 1987; Palla & Zinnecker 1988; 
Ferrara 1998; Uehara & Inutsuka 2002) and the mixing in 
of metals with ionization potentials below 13.6 eV (Dal- 
garno & McCray 1972), thereby initiating the formation 
of large numbers of stars. Similar shock-induced star for- 
mation has been observed and simulated in low-redshift 
intergalactic clouds impacted by radio jets (van Breugel 
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et al. 1985; Fragile et al. 2004). Further, the efficient dis- 
persal of stellar nucleosynthetic products at high redshift 
is required by observations of metals (at levels Z ~ 10~^ 
to 1O~^Z0) in the metal-poor stars of the Galaxy's halo 
(e.g. Freeman & Bland-Hawthorn 2002), in the hot gas 
within galaxy chisters (e..9. Rcnzini 1997; Peterson et al. 
2003), and in the intergalactic clouds producing quasar 
absorption line features {e.g. Tytler et al. 1995; Rauch et 
al. 1996; Songaila 2001). 

Additionally, this formation trigger offers a plausible ex- 
planation of why individual GCs do not have dark mat- 
ter halos today. Previous work (Scannapieco, Ferrara, 
& Broadhurst 2000) has shown that gas in growing pre- 
virialized density perturbations is vulnerable to the influ- 
ence of outflows. Here the dominant mechanism is "bary- 
onic stripping" in which the gas is accelerated above the 
escape velocity and ejected from the associated dark mat- 
ter perturbation. While this has the effect of completely 
suppressing star formation in diffuse pre-viralized regions, 
such interactions may strip the dense gas from virialized 
minihalos while at the same time enhancing their densities 
to form gravitationally-bound clouds. 

In this paper we carry out a two-part study to ex- 
plore this scenario in detail. First, using a simple analyt- 
ical model, we examine heating, cooling, and momentum 
transfer in outflow- minihalo interactions, and determine 
the general properties of star clusters formed by this pro- 
cess. Second, wc use onc-dimensional numerical models to 
examine the turbulence-driven mixing that happens when 
metal-rich supernova ejecta encounters gas of primordial 
composition. Our work builds on recent analyses involving 
one of us (Thacker, Scannapieco, fc Davis 2003; Scanna- 
pieco, Schneider, & Ferrara 2003) in which basic outflow 
properties, global enrichment patterns, and heavy-element 
yields from Population III supernovae were determined. 

The structure of this work is as follows. In §2 we out- 
line a general model for galaxy outflows and high-redshift 
minihalos, and in §3 wc; construct simple analytical es- 
timates for the fate of the minihalo gas subjected to an 
outflow as a function of model parameters. In §4 we use 
one-dimensional numerical turbulence models to examine 
the mixing of metals into the minihalo gas. Conclusions 
are given in §5. 

2. GENERAL FRAMEWORK 

Driven by measurements of the Cosmic Microwave Back- 
ground, the number abundance of galaxy clusters, and 
high-redshift supernova distance estimates (eg. Spergel et 
al. 2003; Eke et al. 1996; Perlmutter et al. 1999) through- 
out this paper we adopt a Cold Dark Matter (CDM) cos- 
mological model with parameters h — 0.7, JIq = 0.3, 
= 0.7, and Qf, = 0.045, where h is the Hubble constant in 
units of 100 km s~^ Mpc~^, and Oq) f^A, and fif, are the 
total matter, vacuum, and baryonic densities in units of 
the critical density, pcrit — 9.2 x 10^'^^ g/cm'^ for our choice 
of h. Note, however, that as the outflow-minihalo inter- 
actions relevant to GC formation are all at high redshifts, 
the value of Oa has no direct impact on our calculations. 

2.1. The Outflow 

To model the expanding outflow we consider a (spheri- 
cally symmetric) Sedov- Taylor blast wave with energy in- 
put cE^Q in units of 10^^ ergs, which is expanding into a 



gas of S times the mean density at a redshift z. While any 
rcaUstic starburst-chiven outflow will show some degree of 
asymmetry, this model gives a resonable approximation to 
the structure of high-redshift starbursts seen in numerical 
simluations {e.g. Mac Low & Ferrara 1999; Mori, Ferrara, 
& Madau 2002; Fujita et al. 2004). Here e is the fraction 
of the total kinetic energy from supernovae (E'ss) which is 
channeled into the galactic outflow. 
In our assumed cosmology, the blast's expansion speed 

is 

(1) 

where 5aa = <5/44, Zs is the redshift at which the shock 
reaches the halo and Rg is the physical (ie. not comov- 
ing) radius of the shock in units of kpc. Our choice of the 
a typical gas overdensity is motivated by the model de- 
scribed in §2.2. The corresponding postshock temperature 
is Ts = 1.4 X 10^ [iis/lOOkm/s)]^ K, assuming an ionized 
gas with a mean molecular weight of = 0.6. By the time 
it reaches Rs, the shock will have entrained a total mass 
of 

M,,total = 1-3 X 106<544 Rl Mq, (2) 

of material and have an overall surface density of 

a, = 1.0 X 10^544 (^^^^ Rs M0kpc-^ (3) 
and it will reach this radius at a time of 

t, = Q.^^5]l,\eE,,)-''-{^-±^Y\l'- Myr. (4) 

We can estimate the total mass in metals as roughly 
2 Mq per 10^^ ergs, a relation that is true both for 
type II supernovae, as well as for pair-production super- 
novae (SN,^-y) from very massive Poplll stars (Woosley & 
Weaver 1995; Heger & Woosley 2002). Assuming that 
half of these metals escape from the host galaxy, we find 
the mass of ejected metals, in units of Mq, to be simply 
Mz = W^E^^. While the Sedov- Tayor solution assumes 
that the entrained material represents swept-up gas, in 
fact a few times the mass in metals is injected into the 
blast wave. However, the Sedov solution will be accurate 
as long as Rg ^ (£'55/5)^/^(1 -|- Zs)~^ kpc, and cooling 
within the bubble is small. 

To estimate the energy in a typical high-redshift out- 
flow, consider a young galaxy whose total mass is Mgai > 
10^ Mq. Such a galaxy will have a virial temperature ex- 
ceeding 20,000 K, so its gas readily cools via atomic (HI) 
line emission, leading to star formation. The assumption 
that 10% of the gas is converted to stars gives good agree- 
ment with observed high-redshift star- formation rates, as 
well as with abundances of metals measured in high- 
redshift quasar absorption line systems (Thacker, Scan- 
napieco, & Davis 2002; Scannapieco, Ferrara, & Madau 
2002). For very massive (PopIII) stars, the models of 
Heger and Woosley (2002) suggest that supernovae pro- 
duce some 10^^ ergs for every 30 Mq of material in new 
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stars. For less massive (PopII) stars, a Salpctcr IMF yields 
one supernova with an explosion energy of 10"''^ ergs for 
every 150 Mq in new stars (eg. Tegmark, Silk, & Evrard 
1993). Thus we obtain 

i;5r,(PopIII) ~ 50Af9 and £5r,(PopII) ~ lOMg 

(5) 

where here and below Mg is the galaxy mass in units of 
10^ Mq. For the wind efSciency (e) we rely on the simula- 
tions described in Mori, Ferrara, & Madau (2002), which 
indicated e ~ 0.3 in the case of a 2 x IO^Mq star-bursting 
galaxy. 

Finally, we note that a lower bound, eE^^ > 1, follows 
from the fact that high-redshift outflows dispersed metal 
efficiently enough to preclude stars of primordial composi- 
tion being formed today (Scannapieco, Schneider, & Fer- 
rara 2003). This minimum value, plus the numbers cited 
above, suggest a fiducial outflow model in which e = 0.3 
and i?55 = 10. 

2.2. The Protocluster 

Having established basic outflow parameters, we now 
develop a simple model for the gas and dark matter in a 
protocluster whose total mass = Mq x IQ^Mq. We 
assume that the gas has a primordial composition (76% H 
and 24 % He, by mass) and is unionized, giving a mean 
molecular weight /i = 1.2. 

At a redshift Zc ~ 10, corresponding to a cosmic age 
1/2 Gyr, the gas and dark matter collapse and virial- 
ize. Initially, the mean density of the protocluster cloud is 
enhanced by a factor A above the background, 



Pc = AOo(l + ZcfpcTit, 



(6) 



where the enhancement factor for a virialized cloud at high 

redshift is well approximated by the value in a critical uni- 
verse A = 178 {e.g. Eke, Navarro, & Frenk 1998). With 
this choice, the cloud's virial radius is 



1/3 



Rr = 0.3 M, 



and its virial velocity is 



1 



10 



kpc. 



(7) 



Vc = 4:.4:M' 



1/3 / 1 + Zc 



10 



1/2 



kms-^ (8) 



As first observed in the numerical simulations of 
Navarro, Frenk, and White (1997, hereafter NFW) we as- 
sume that the CDM minihalo develops a radial profile of 
the form 



p{R) 



floPc 



cx{l + cxf 3F{c) 



(9) 



where x = R/Rc, c is the halo concentration parameter, 
and 

i^(i) = ln(l+t)-^. (10) 

We then assume that as the procluster gas collapses within 
the dark mater halo, it is shock heated to the virial tem- 
perature 

Te = 720M6'/3l±^ K, (11) 



and develops the density distribution of isothermal matter 
in the CDM potential well: 



(12) 



The central gas density (po) is determined by the condition 
that the average baryonic density within the virial radius 

is equal to (ilb/^o)Pc- The escape velocity for an atom at 
a radius R within the well is given by 

vUR = xRc^) = 2f^M2^dR' 



■ 2v: 



F{cx) 



1+cx 



xF{c) 



(13) 



so that vl^M = 2vlc/F{c). 

Following Madau, Ferrara, & Rees (2001) we assume a 
typical concentration parameter of c = 4.8, although there 
are some indications that high-redshift halos may be less 
concentrated than expected from this estimate {e.g. Bul- 
lock et al. 2001). With this choice p{Rc) = 4406Pcrit, 
motivating our typical value of ^44, and Wesc(O) = 10 ^c> 



%sci-^c) — 3.9 w;:, and 



Po 



r^{i+t)^iHut 



Pcrit = 1970pcrit, 



(14) 



where t = cx bxvA A = 2c/F{c) = 10.3. 

Note that a gas-to-star conversion efficiency of 10% fol- 
lowed by a loss of 50% of a young globular cluster's orig- 
inal stellar mass (through winds, tidal forces, etc.), when 
combined with the aforementioned present GC mass limit, 
< 1O^M0, implies that the likely limiting mass of 
proto-globular clouds (gas-l-GDM) is 



Mg (protocluster) < 10^ 



(15) 



For a collapse redshift Zc ~ 10, the corresponding limit 
on the virial temperature is Tc < 15, 000 K. Choosing a 
20% conversion efficiency gives a corresponding limit on 
the virial temperature of Tc < 9, 000 K, which is low- 
enough to preclude substantial ionization, hydrogen line 
cooling, and "unassisted" star formation. Thus, given the 
uncertainties involved, the maximum size of minihalos cor- 
responds to the maximum size of halo globular clusters. 

Lately, a great deal of attention has been given to the 
inner profile of dark-matter halos, as measured from the 
properties of low surface brightness (LSB) galaxies, whose 
gravitational potentials are thought to be dark-matter 
dominated at all radii. In particular, studies of these 
objects have found that their inner rotation curves are 
likely to be well-described by constant density cores, in 
contrast to the \/R profile assumed here (e.g. Carignan 
& Beaulieu 1989; Carignan, & Sancisi 1991; de Blok & 
McGaugh 1997). Recent studies have complicated the is- 
sue, however, suggesting that to some degree, these re- 
sults may have been affected by the poor angular reso- 
lution of the H I observations. By accounting for beam 
smearing, several groups have shown that the H I rotation 
curves of LSB galaxies are consistent with a wide variety 
of dark matter potentials, ranging from constant density 
cores to profiles as steep as R~^ {e.g. Blais-Ouellette et 
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al. 1999; Swaters, Madorc, & TrcwhcUa 2000; van den 
Bosch & Swaters 2001). Further high-resolution Ha mea- 
surements have reached similar conclusions (Swaters et al. 
2003). While still an unsettled issue, the presence of such a 
core will only reduce the gravitational potential in the cen- 
ter of the cloud, enhancing the impact of a cloud-outflow 
interaction. This effect is only important in the interior 
R < Rc/c profile of the cloud, however, which contains 
~ 15% of the gas mass. Thus we expect our NFW model 
to provide a reasonable description of the the protoclustcr; 
at worst it provides a slight underestimate of the impact 
of shocking by galaxy outflows. 

Finally we note that, given our cosmological model (with 
a primordial power spectrum as given by Eisenstein & Hu 
1999, (Jg = 0.87) and the efficiencies assumed above, by 
Zc = 15 the region that evolved into the Milky Way was 
filled with ^ 1000 minihalos that could form into GCs with 
stellar masses > lO^M©; by Zc = 10, there were ~ 3000 
such minihalos. These estimates are based on the Lacey & 
Cole (1993) progenitor model and an assumed total Galac- 
tic mass and formation redshift of 2 x IO^^Mq and Zf = 2. 
Even assuming just the mean cosmological density results 
in ^ 500 and ^ 2000 such minihalos at Zc = 15 and Zc = 10 
respectively (Press & Schechter 1974). Thus, although we 
show below that only a fraction of minihalos will be located 
the correct distance from a starburst to form a GC, our 
mechanism need not be particularly efficient to account 
for the Galaxy's 100 halo globular clusters. 

3. FATE OF THE PROTOCLUSTER 

3.1. Outflow- Protoclustcr Interactions 

Using the simple outfiow and protoclustcr models de- 
scribed above, we now outline the general features to be 
expected when the shock interacts with the gas in a CDM 
minihalo. The mass of the shell that impacts such a pro- 
toclustcr is 

Ms = ttRIcts = 2.9 X 10* (544 M^^^'^w^^^^Rs Mq, 

(16) 



10(l + Zc)2 

and at the time of impact its momentum is 



Ps = ttRIps = 2.2 X 10^ Ml'"" ((544 e ^55)'/^ 

>< '°7iK/'^ -^«"'^'^0kms-\(17) 

where ps is a momentum surface density. Recalling that 
Mz = IQ'^E^^Mq is the total mass in metals within the 
blast-wave material, it is straightforward to calculate that 
the metallicity to which the protoclustcr gas is enriched is 

-2 



= 0.076 ^Mg ^'^Ey 



55 



10 



R„ 



Zq, (18) 



where we have taken solar metallicity to be 1/50 by mass, 
and have assumed that a fraction ^ of the metals is thor- 
oughly mixed into protoclustcr gas and is contained in the 
GC stars observed today. This mixing eflaciency is studied 
in detail in §4. 

We are most interested in clouds whose stars are en- 
riched to approximately 1/30 the solar value (Ashman & 
Zepf 1998), which occurs at a typical distance of 

-1 



Rs = 1.5 M, 



-1/6 



(Ci^55) 



1/2 



10 



kpc. (19) 



This is sufficiently large that the original mass input is 
negligible compared to the entrained mass (see eq. 2). The 
time it take for the bubble to reach this radius is 



ts = 1.5 e-1/2 ^5/4 ^3/4 



44 

x(^r(^)" 



■5/2 



^r/so' Myr. 



(20) 



Typically, this is shorter the cooling time within the bub- 
ble 



^cool — 



= 18 Te^. 



44 



1 f'^ + Z 

To~ 



A:^3 Myr, (21) 



where A_23 is the radiative cooling rate of the gas in units 



of 10- 



ergs cm s and is electron number density 



of the shocked gas in units of cm~^, Tq is the gas temper- 
ature in units of 10^ K, and we have asssumed that the 
density of the post-shock gas is enhanced by a factor of 4. 

Thus our assumed Sedov solution should be reasonable for 
the range of values considered in this study. 

Futhermore, for likely values of 2^55, Rg is compara- 
ble to the virial radius of the outflowing galaxy, which is 

30Mg^'^(l -|- Zs)^^ kpc. Hence, the density of the medium 
through which the shock is passing is substantially higher 
than the mean IGM density, and more like our fiducial 
value, i5 = 44, taken at the virial radius. 

From Eqs. (1), (2), (17), and (19) we determine the 
shock velocity at the time of impact to be 



1 +Zc 
1 + 2:. 



3/2 



kms ^, 



(22) 

and the mass and momentum impinging on the cloud to 
be 



and 



4.3xlO'{6lMe^E,,y/'(^^±^y Mq, (23) 



= 1.8 X 10^ {S^Me £2 ^-1 £55)1/4 



1 + Zs 
l+Zc 



3/2 



Mq kms-^ (24) 



If we compare Pg with the virial velocity of the cloud times 
its baryonic mass. 



Pr = 6.3 X 10^M« 



4/3 / 1 + 



10 



1/2 



Mq kms-\ (25) 



we see that for redshifts ~ ~ 10 the blast wave's mo- 
mentum is sufficient to move the gas out of the dark matter 

potential well whenever Mg < 10 ( — 1 . But, how 

will such a dark cloud be impacted as a function of radius, 
and what is the ultimate fate of its swept-up gas? 
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3.2. Three Stages of Evolution 

In order to study this interaction in more detail, we 

adopt a specific set of parameters for §3.2, and in §3.3 wc 
discuss the effect of varying these values. Here we put E^^, 
e = 0.3, Zc = 10, Zs — 8, and Zc = 1/30^0. We consider a 
protocluster with M = 3.2 x lO^M© whose total gas mass 
is approximately 5 x IO^Mq. Assuming an overall star 
formation efficiency ^ 10%, this provides a good match 
to the observed peak in the globular cluster mass function 
at IO^Mq (in stars). In this section, we assume that all 
metals eventually mix into the protocluster gas, so ^ = 1. 
Altogether, these choices represent our fiducial model. 

We identify three important stages in the evolution of 
shocked cloud. The first of these occurs as the outflows 
moves across the minihalo. At this point, the key question 
is whether the impinging momentum is sufficient to accel- 
erate the gas to its escape velocity, stripping it from the 
dark matter potential. To estimate when this occurs, we 
compute the average momentum surface density as 



Pcib) = 2 



de 



(26) 

b is the impact parameter from the central axis of the 
cloud, and £ is the distance along a line parallel to this 
axis. This is plotted in third row of Fig. 1, for our fiducial 
model. 

Comparing the momentum surface density with the 
momentum per unit area in the outflow, we find that 
Ps > Pcib) for all impact parameters ^strip — 0.2i?c. 
In our simple picture, we assume that all the gas out- 
side of 6strip will be stripped from the potential, while the 
denser central regions will be left behind, resulting in an 
elongated bell-shaped distribution, not unlike the coma of 
a comet. 

In order to better relate our radial profiles to the prop- 
erties of the cloud, in the upper two panels of Fig. 1 we 
plot 

M{> 6) =277 / db'adb'), (27) 



the mass outside a given impact parameter, and 



Zcib) 



ZcM 
nRl<Jc{b)' 



(28) 



the local metallicity at that distance, where adb) = 

2 ^ dzpgg^siVb^ + z"^). From these estimates, we 

find that approximately 80% of the gas mass is able to 
be efficiently stripped from the minihalo and that, once 
homogenized, the metallicity of this gas exceeeds IQ^^Zq. 

The second important stage occurs just after the blast 
front passes across the protocluster. At this time, the gas 
has been heated to a postshock temperature of several mil- 
lion degrees and its density has been enhanced by a factor 
of four. There are three time scales that then enter into 
the problem. The first of these is the sound crossing time, 
which we estimate as 



_ {Rl-b^Y/^ 
4c., 



(29) 



1/2 

where Cg = OASvg = 115Tg km/s is the postshock sound 
speed. 

The first key issue is then whether the self-gravity of 
the gas alone can resist the pressure associated with such 
an enormous temperature increase. This is determined by 
comparing the sound crossing time with the free-fall time, 
which can be written in convenient units as 



ts 



32Gp, 



67 n' 



-1/2 



Myr, 



(30) 



gas 



:6.8r-i/X^'(i + ^c)-' [1 



{b/R,ff^ Myr, 



where G is the gravitational constant, and n is the den- 
sity of the gas in atoms cm^''. In the fourth panel of Fig. 
1 we see that for our fiducial set of parameters tg ^ tsc 
at all impact parameters. Thus it seems the swept away 
gas will be evaporated into the IGM within a sound cross- 
ing time. In order for this to occur, however, the gas 
must expand before cooling process arc able to dissipate 
its thermal energy. This is determined by the cooling time 
as given by eq. (21). Note that in this equation A_23 is a 
function not only of the temperature of the gas, but also 
of its elemental composition, as the presence of heavy ele- 
ments greatly increases the number of transitions that can 
radiate efficiently. For simplicity we estimate this radia- 
tion by assuming that all the gas is at the mean metal- 
licity given by cq. (18) and by taking solar abundance 
ratios, which allows us to make use of the tabulated mod- 
els of Sutherland & Dopita (1993). This approximation 
is dependent on the prompt mixing of the minihalo gas 
with all the impinging material (as discussed in §4), but 
in fact a range of mctallicities may be found in the gas 
at a given impact parameter. For temperatures within 
lO^K > T ^ 10"^, this leads to A-23{Z = 1/30Zq) values 
ranging from 0.6 to 10. At higher temperatures, cooling 
is dominated by bremsstrahlung and is largely metallicity 
independent, while molecular cooling becomes important 
at lower temperatures, as discussed below. Finally we es- 
timate the post-shock density at an impact parameter b as 

a^b) [(i?2- 62)1/2/4] 

With these simplifying assumptions, wc obtain the cool- 
ing times that are plotted as the dashed lines in the fourth 
panel of Fig. 1. Note that fcooi is several orders of magni- 
tude smaller than the free-fall time, as can be inferred di- 
rectly from the properties of GCs observed today (Murray 
& Lin 1992). In fact, for a large range of impact param- 
eters, fcooi is much smaller than the sound crossing time. 
Thus, despite enormous postshock temperatures, efficient 
radiation by the dense and metal-c;iiric;hed halo gas is able 
prevent evaporation of the gas within b < 0.7 Rc, which 
we label as 6cooi- This means that roughly half of the gas 
(0.2i?c ^ b < 0.7i?c) is expelled from the dark matter 
potential, yet cools sufficiently quickly to remain gravita- 
tionally bound. Modulo our assumption of efficient mixing 
(^ ~ 1), this gas has a relatively small range of metallici- 
ties, from about 10~^ to 10~^-^ Zq. 

At high temperatures, the cooling time is a monotonic 
function that strongly decreases with decreasing tempera- 
ture. The cloud becomes ever more efficient at radiating 
its energy until the gas begins to cool below lO^K, when it 
reaches the third and final important stage. At this point. 
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the gas is largely neutral, and the cooing rate decreases 
precipitously. Two new processes then become important, 
the first of these is the production of molecular hydrogen 
by nonequilibrium reactions {e.g. Ferrara 1998) and the 
second is infrared line emission by C ii, Fe ii, and Si ii, 
whose ionization potentials are less than 13.6 eV (Dalgarno 
& McCray 1972). In the case of molecular hydrogen, run- 
away cooling from ~ lO^K results in appreciable levels of 
H~ and H^, which act as intermediaries in H2 formation. 
While the exact numbers are uncertain to within a factor 
~ 3 this can result in levels of H2 ^ 1% (see e.g. Uehara 
& Inutsuka 2000, Figure 1), which, in turn yield coohng 
rates of ~ lO-^^T^-^ ergs cm^ s"! (Galli & Palla 1998). 
The atomic infrared cooling is somewhat less efficient but 
also a weaker function of temperature, ranging from about 
10-27 ergs cm3 s'^ at 1000 K to lO^^s ergs cm^ s'^ at 100 
K for 1/30^0 gas. Thus cooling through H2 is expected 
to be slightly more important than metal line cooling at 
1000 K and slightly less important at 100 K. 

The relevant times scales at late times arc plotted in the 
bottom panel of Fig. 1, assuming that the gas cools at a 
fixed density. Here we see that at both 1000 and lOOK, 
the cooling time is substantially shorter than the sound 
crossing time at all radii, and the gas is likely to cool to 
very low temperatures even if photo-dissociating radiation 
(which we do not attept to model) were strong enough to 
quickly destroy the formed H2, meaning that all cooling 
was through C 11, Fc 11, and Si 11, or if mixing is ineffi- 
cent, meaning that all cooling was through H2. In fact, at 
100 K the sound crossing times even lie above the free-fall 
times. Thus, although the cooling rate is orders of mag- 
nitude smaller at these low temperatures, weak radiation 
from H2 or metals is still sufficient to cool the cloud to 
the point at which its thermal pressure cannot counteract 
self-gravity. Note that this runaway collapse takes place 
even if we assume no additional density enhancement dur- 
ing cooling, and no external pressure from the hot medium 
that is likely to be found behind the shell [as seen for ex- 
ample in NGC 3077 (Ott, Martin, & Walter 2003)]. As 
ts oc and tcooi oc n~^, while tgc oc R (x n~^^^, such 

density changes during cooling from ^ lO^K will only en- 
hance collapse at low temperatures. Furthermore external 
pressure, while negligible for the hot cloud, may equal or 
even exceed the thermal pressure of ~ lOOOK gas. Thus 
we believe low-temperature cooling does not represent an 
important barrier for star formation in shocked minihalos. 

To summarize, the two key issues that determine the 
fate of the gas are momentum transfer to the cloud dur- 
ing outflow shocking (stage 1) and the ratio of the cooling 
time to the sound crossing time just after the shock moves 
across the cloud (stage 2). It is these issues that we now 
address in some detail, as we study the impact of varying 
several model parameters about their fiducial values. 

3.3. Dependence on Input Conditions 

Having outlined the general stages we expect shocked 
minihalos to undergo, we now consider a range of values 
for the input parameters of our model, and relate these 
to the observed range of properties of halo globular clus- 
ters. The first issue we explore is the impact of varying the 
distance between the minihalo and the outflowing galaxy. 
As this is directly related to the mass in metals reaching 
the cloud, according to eq. (18), this constrains the range 



of globular cluster metallicities that can be generated in 
our model. Holding all other fiducial paramters fixed, we 
select mean values of IQ-^ o, IQ-^-^, lO'^, IQ-^-^, and 
10-^^0, which correspond to distances Rs of 21, 12, 6.6, 
3.7, and 2.1 kpc, respectively, and to the results that are 
plotted in Fig. 2. 

At the smallest distance, more momcintum is imparted 
to the halo, causing stripping to be slightly more efficient. 

— 1/2 

As the momentum surface density is proportional Rs 
however, this effect is minor in comparison to the change 
in the postshock temperature Ts oc oc Rj^. Thus, in 
this case, tcooi is vastly increased while tsc is decreased, 
and the stripped gas is evaporated without forming stars. 
Similarly, increasing the distance has a much stronger im- 
pact on tsc and tcooi than on ps . In the models in which 
Z = 10-2, 10-2-5, and 10~^ Zq, almost all of the gas out- 
side of 6 w 0.3, 0.35 and 0.4i?c respectively is stripped 
from the halo, yet able to efficiently cool, which again cor- 
responds to roughly 50% of the gas mass. In fact the max- 
imum distance at which ejection and shock-induced star 
formation are effective is not determined by the stripping 
criteria, but rather by the postshock temperature being 
too low to ionize the minihalo gas and allow it to radiate 
efficiently. The post-shock temperature for Zc = 10-^.^0, 
for example, is only ~ 5000K, well below the lO^K 
needed for effective collisional excitation of hydrogen. 

These results compare well with the observed metallicity 
distribution of halo globular clusters, which peaks slightly 
below [Fe/H] = —1.5 and falls off gradually with very few 
objects having iron abundances below [Fe/H] = —2.5 or 
above [Fe/H] = —1 (Ashman & Zepf 1998). Furthermore 
the relatively constant metallicity between 0.2Rc <, b <j 
0.7Rc compares well with the observed A[Fe/H] < 0.1 dex 
spread observed within individual objects {e.g. Suntzeff 
1993). 

The second important issue is the dependence of our 
model on the mass of the shocked cloud. In Fig. 3 we vary 
this quantity over two orders of magnitude, while keeping 
other fiducial parameters fixed. By eq. (18), this means 
that the distance between the source galaxy and the mini- 
halo is assumed to increase with decreasing mass. Thus, 
although stripping is seen to be more efficient in the lower 
mass case, the post-shock temperature is reduced too, as 
Tg (X Mc^^ at a fixed Z,,. This greatly reduces the 

cooling time, allowing the gas as a whole to be stripped 
and to cool quickly enough to remain bound. 

From this comparison it is clear that our assumed mech- 
anism becomes more effective in the case of smaller mini- 
halos, and that the minimum scale of globular clusters 
probably depends instead on processes that occur after 
star formation. Fortunately, as discussed in §1, this limit 
is easily understood in terms of the dominant destruction 
processes that act on GCs in the Milky Way, in particular 
disk-shocking and mechanical evaporation. On the other 
hand, the maximum mass of GCs is likely to be a feature 
of the proto-cluster cloud itself. Turning our attention to 
the high mass case, we find that if Mc = IO^^'^Mq then 
^'strip = 0.5iic and 6cooi = 0.25i?c. Thus we expect our 
scenario to be relatively inefficient in forming stars in ob- 
jects more than a few times larger than our fiducial mass, 
Mc = W^'^Mq, which rougly corresponds to the largest 
minihalo as well as the maximum observed GC mass, as 
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discussed in §2.2. 

Next wo evaluate the effect of varying the energy input 
parameters £^55 and e. Increasing £^55 increases the mass of 
ejected metals, so raising this value also raises the distance 
corresponding to 10~^'^Zq enrichment. Thus, changing 
the number of supcrnovae in the outflowing galaxy has a 
relatively minor impact on the momentum imparted to the 
cloud, changing fostrip by only a few percent, as is shown in 
the left panels of Fig. 4. Similarly, simultaneously increas- 
ing the energy input and distance serves to damp the effect 
has on the sound crossing time and and cooling time of 
the postshock gas. Note that as Rg oc E^/^ the postshock 
temperature is proportional to Rj^E oc decreas- 
ing slightly with increasing £^55, and shorting icooi while 
raising t^c somewhat. Thus the primary impact of £^55 on 
our globular cluster formation model is simply to shift the 
relevant distance further away from the outflowing galaxy, 
while keeping the metallicity fixed, and slightly improving 
the fraction of minihalo gas available for star formation. 

In the right columns of Fig. 4 we study the impact of 
varying e, the ejection efficiency. As this parameter affects 
only the shock velocity, it has no impact or our choice of 
fiducial distance and is not subject to the damping ob- 
served in the left panels. In this case 6strip ranges from 0.2 
Rc to 0.35 Rc while 6cooi ranges from 0.4i?c in the e = 0.5 
case to iic if e = -1. Thus the wind efficiency represents a 
larger model uncertainty than £^55 , although at least 40% 
of the gas mass lies between fegtrip and 6cooi in all cases. 

Finally, in Figure 5 we study the effect of varying the 
density of the medium through which the shock is mov- 
ing, as determined by S, and the density of the minihalo, 
as determined by Zc- From eqs. (22) and (24) ps oc S^^^ 
while Tg oc S^^. Thus increasing S improves the efficiency 
of our formation mechanism, pushing 6strip towards the 
center, while shifting bcooi outwards. Similarly, decreas- 
ing S reduces both the efficiency of stripping and cool- 
ing. At i5 = 10, 6cooi = 0.15i?c < 6strip = 0.35i?c and 
our GC formation mechanism fails completely. Thus, our 
scenario is most likely to take place at separations that 
are comparable to the virial radius of the starbursting 
galaxy. This point is particularly noteworthy because, 
while protocluster-scale minihalos are ubiquitous in CDM 
cosmological simulations, no patently intergalactic globu- 
lar clusters have been found to date (McLaughlin 1999; 
Mari'n-Franch & Aparicio 2003). 

Changes in the density of the cloud, as determined by its 
formation redshift, have the opposite effect on the shock 



properties. In this case Ps cc {I + Zc)' 



-3/2 



-1/2 



and 



Tg oc (1 + Zc)^ oc Sc- Variations in Zc also effect the prop- 
erties of the cloud, however, in particular shortening the 
free-fall time and increasing the imparted momentum Pc 
as the formation redshift is increased. Thus later- forming 
clouds are heavily favored by our mechanism, and is un- 
likely that efficient star-formation was achieved in objects 
in which Zc^ Zg. Note that it is the ratio of 1+Zs and l+Zc 
that is important for this comparison, and thus we expect 
little difference in our results if both the collapse and out- 
fiow rcdshifts were shifted to higher values, as might be 
necessary if cosmic reionization took place early, as sug- 
gested by polarization measurements from the Wilkinson 
Microwave Anisotropy Probe (Kogut et al. 2003). 

To summarize these sensitivity studies in the context 



of our scenario, cloud-outflow interactions are efficient at 
forming stars in clusters with mean metallicities rang- 
ing from about IO'^-^Zq to 10'^ Zq. Below IO'^-^Zq 
shocks are too weak to induce star formation, while above 
1O~^Z0, shocks are too powerful and disrupt the gas com- 
pletely. Induced star formation is largely independent of 
the total energy input from a given starburst, although 
some uncertainty is introduced by the unknown fraction 
of this energy that goes into powering a galaxy outflow. 
While outflowing shells arc efficient at forming stars in 
small minihalos, this become more difficult in the largest of 
such objects, whose gas masses compare well with the max- 
imum stellar masses of globular clusters. And finally, our 
mechanism is most efficient in cases in which the density 
contrast between the cloud and the surrounding medium 
is not too high: in minihalos close to starbursting galaxies 
that formed at a similar redshift. 

4. TURBULENCE MODELS OF METAL MIXING 
4.1. A ID Code for Turbulence Transport and Mixing 

In §3 we saw that the chemical homogeneity observed 
within globular clusters can be associated with the range 
of metallicities generated in outflow-minihalo interactions. 
Yet these estimates were based only on the mass of heavy 
elements impinging on the cloud as a function of impact 
parameter. For true homogeneity within the resulting stel- 
lar population to be achieved, not only must sufficient 
masses of metal be present, but this material must be well 
mixed into the star- forming gas. Note that this chemical 
structure is of particular importance as it plays a key role 
in gas cooling as described in §3.2. 

As a first step towards examining this "fine grained" 
mixing (Rees 2003) in detail, we have developed a one- 
dimensional numerical code, which make use of a well- 
tested generalization (Besnard et al. 1992) of the the 
widely used "iC— e" model of turbulence transport (Harlow 
& Nakayama 1967). In this case the equations of conser- 
vation of mass, momentum, and energy are coupled to the 
Reynolds stress tensor, Rij = pu'-u'- which represents the 
ensemble average of the product of the density and the 
velocity departures from the mass-weighted mean. Thus 
we have 



dp _^ dpui 
dt dxi 
dpUi dpuiUi 



= 0, 



at 

dpE 



+ 



+ 



dxi 
dpEui 
dxi 



dR^^i 
dxi 
dRi 



dP 

dxi 

dPui _ 
dx\ dx\ 
d f pvtCp dT 

dxi V Pr* dxi 



(31) 
(32) 

(33) 



where t and xi arc time and position variables, p{x\,t) 
is the mass density field, Ui{xi,t) is the mass-averaged 
mean-flow velocity field in the i direction, P{xi,t) is the 
mean pressure, E{xi,t) is the total mean energy, T{xi,t) 
is the mean temperature, Ui{xi,t) is the volume averaged 
mean-flow velocity field in the i direction, and we take 
Prt = 1.0 and Cp = 3/2, as appropriate for a polytropic 
gas. The Reynolds stress tensor is then modeled using a 
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Boussinesq-type approximation as 



dxj 



+ 



duj_ 
dxi 



2p 
3 



duk 
dxk 



+ K], (34) 



where Vt is the turbulence "eddy" viscosity (which is calcu- 
lated as Vt = K{xi,t) is the turbulence energy 
density, and e{xi,t) is the turbulence energy density decay 
rate. The transport equation for each of these quantities 
is constructed as a compressible generalization of the stan- 
dard K — e model, namely 

dpK dpKui _ ^ dP ^ dui ^ d fdK 
dt dx\ dxi ^'^ dxi dxi \^''dxi 

(35) 

where a, = p'u^/p (with p' and u' the departures from 
the volume- weighted mean density and velocity), and the 
terms on the right hand side of the equation represent 
buoyancy creation, shear creation, decay, and self-diffusion 
respectively. Similarly 

dpe dpeui e \ dP du, ' 

at oxi K I oxi axi 

dxi \ (Te dxi J ^ dxi ' 

where the first four terms on the right hand side parallel 
those in eq. (35), and the final term represents changes in 
the turbulent scale due to expansion or contraction of the 
fluid element. This scale can be defined as 



^3/2 



(37) 



such that an expanding motion results in an increase in S, 
which translates to a decrease in e, likewise a contraction 
decreases S, driving e up. In eq. (36) pai is the turbulence 
mass flux, while the empirically fit constants are 1.3, 
Ce, = 1.55, Ce, = 2.0, Ce, = i + and C,, = 1.25, 

{e.g. Launder, Reece, & Rodi 1975) where d is the dimen- 
sionality of the compression modeled in the last term of 
eq. (36), which in our case is 1. Closure of our system of 
equations is achieved by the transport equation for af. 



dpai dpaiUi dp 



dt 



■ + 



dxi 



p dxi ' 



pe btrh dP 

-Cai—ai-\ -— , (38) 

K p oxi 



where the first two right-hand terms represent turbulence 
distributive creation and "drag" decay, C'ai = 2.0, and the 
last term with 6trb = p'p' /p- describes buoyancy effects, 
which are here neglected. Mixing between multiple mate- 
rials is computed by carrying additional transport equa- 
tions for their concentrations, c^. For this purpose we use 
a diffusional approximation for the species flux: 



dpCa dpCaUi 



dt 



dxi 



d_ 

dxi 



where a is an index over materials and ctc = 1.3. We 
compute Ui as Ui — at and the temperature as T = 

^3^" (^E — K — ij-^ , with ks the Boltzmann constant. 

Our code does not attempt to account for self-gravity or 
the impact of the dark matter gravitational potential, as 



these only become of primary importance when the full 
two-dimensional geometry of the distortion is calculated, 
as discussed in §4.2. 

This formulation of turbulence equations is closely re- 
lated to that used in the CAVEAT numerical code (Adessio 
et al. 1990), where it has been compared with a wide range 
of experimental results. Our implementation is an Eu- 
lerian finite-difference one, with donor-cell advection be- 
ing used throughout. Test calculations show overall qual- 
itative agreement with experimental results for the shear 
problems of interest us. While our approach is sufficiently 
accurate to demonstrate our ideas for this paper, more 
precise analyses may require the use of a full spectral rep- 
resentation, like that described by Steinkamp, Clark, & 
Harlow (1999a; 1999b). 

4.2. Application to Shocked Minihalos 

Having developed and tested a simple one-dimensional 
code to track turbulent mixing in compressible fiuids, we 
then applied this tool to the problem of cloud-outflow 
interactions. Here we focused on two key geometries: 
an edge-on configuration, in which mixing occurs due to 
Kelvin-Helmholtz instabilities as the outflow shears along 
the side of the cloud, and a face-on configuration, in which 
mixing due to Raylcigh- Taylor instabilities occurs as the 
outfiow material moves into the denser cloud. Each of 
these configurations is illustrated in Figure 6. 

Turning our attention first to the shear instability case, 
we chose our initial conditions to approximate the fiducial 
problem described in §3. Our simulation domain was di- 
vided into 800 evenly-spaced zones, spanning a distance 
of 840 physical parsecs, twice the radius of our fiducial 
M = 10^-^Mq,Zc = 10 cloud. The first 400 of our simula- 
tions zones were filled with material 1 , which was taken to 
represent the primordial gas making up the cloud. Here 
the density was computed as per eq. (12), the fiuid was 
initially at rest, and the temperature was taken to be the 
virial value of 1650 K. The second 400 zones were filled 
with material 2, which was taken to represent the enriched 
outfiow. In these cells the density was taken to be 4 times 
the density at the radius of the cloud, Uy was 200 km/s 
(~ 200 parsecs/Myr), and the temperature was taken to 
be the postshock value of 560, 000 K. On both sides of the 
shock, material is assumed to be largely quiescent, with 
an initial turbulent kinetic energy equal to only 2% of the 
total energy. An initial turbulent scale S of 10 parsecs was 
chosen in both materials. 

These conditions are shown in the various panels in the 
left column of Figure 7. In order to better quantiiy mixing 
in the simulation, we also computed 



M{b) = 1 



Jq dxci {x)p{x)x'^ 
J^'' dx"ci{x')p{x')x''^'' 



(40) 



a quantity that estimates the fraction of the total cloud 
gas mass contained outside of a given impact parameter b, 
which is plotted as the dashed lines in the upper panel. 

Adopting a time step of 600 years and a time unit such 
that f = 1 is defined as Rc/vs = 1.5 million years, we ran 
our simulation to a final time of f = 5. From the estimates 
in §3.1, one can expect a ^ 20 Myr delay between shock- 
ing and star formation, which is primarily governed by the 
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free-fall time; however by these late times the distortion of 
the cloud is likely to be so severe that our one-dimensional 
calculations will look nothing like the true configuration. 
Thus we choose to stop our simulation much sooner, quot- 
ing our simulated mixing as a lower limit. These results 
are depicted in Figures 7 and 8. 

Three major features are visible in the mean-flow quan- 
tities represented in the upper panels in these plots: a 
strong shock moving into the cloud, a material interface 
moving inwards, and a rarefaction wave moving outward 
into the enriched gas. The motion of each of these fea- 
tures is well approximated by the analytical solution for 
the propagation of a strong pressure discontinuity in the 
absence of turbulence transport or shear (eg. Harlow & 
Amsden 1971): 

'^x,int 3(^s, outflow *^5,r)? (^-^) 

4 

Ux,s ~ -^Ux,int, and (42) 

'^x^r ^ ^5, outflow; (^*^) 

where Wx.int, Ux,s, and Ux,r are the velocities of the inter- 
face, shock, and rarefaction wave respectively, Cg^outflow is 
the sound speed in the outflow region, and Cs^r is the sound 
speed to the left of the rarefaction wave. Initially, the in- 
ward shock and fluid interface move at speeds of ~ 90 and 
~ 65 km/s respectively, but slow somewhat as they reach 
the dense regions near the center of the cloud. On the 
other hand, the sound speed in the unshocked regions is 
quite small (^ 5 km/s). As the cloud is assumed to be 
initially in pressure equilibrium, any motions due to grav- 
itational effects should also be ~ 5 km/s, justifying our 
neglecting self-gravity and the dark-matter gravitational 
potential in this one-dimensional calculation. 

Turning our attention to the turbulence energy profile, 
we find a notable increase along the shearing interface be- 
tween the two materials. This is initially confined to a 
relatively narrow region at early times, but expands to 
cover ~ 100 parsecs hy t = 5, and the corresponding tur- 
bulent diffusion is sufficient to mix metals into ^ ~ 20% of 
the cloud by this time. In order to distinguish this mix- 
ing from the numerical diffusion intrinsic to our Eulerian 
approach, we also include values in the right column of 
Figure 8 for an identical run in which Vt has been set to 
zero in eq. (39). 

As a test of convergence, both runs were repeated dou- 
bling the resolution and halving the overall time step. This 
had the effect of reducing the mixing in the vt = Q esti- 
mates of the concentrations, while leaving other quantities 
unchanged in both runs. 

To explore our model in the face-on case, we conducted 
a second simulation with initial conditions as shown in the 
left column of Figure 9. Here the simulation domain was 
divided into three regions, the left two representing the full 
density profile through the cloud, and the right represent- 
ing the impinging outflow. The simulation domain was 
divided into 1200 cells, each with the same width that was 
used in the edge-on calculation (1.05 parsec). Again the 
dashed lines in the upper panels represent the total mass 
of primordial material contained in front of a given i value, 
M{t), where i is defined as the distance from the back of 
the cloud. This mass is calculated using an equation anal- 
ogous to eq. (40), but now adding over the full cloud such 



that M ~ 0.5 at the center of the density profile. 

Choosing a time step and time unit as in the edge-on 
case, we ran our simulation to a final time of f = 5, with 
results as depicted in Figures 9 and 10. Again, three major 
features are visible in the mean-flow quantities. As in the 
edge-on case, a strong shock and material interface quickly 
propagate inwards. Initially, the fluid interface moves at 
the incoming speed of 200 km/s, and Ux,s and Ux,int share 
roughly the relation expected from eq. (42). But as in 
the Kelvin-Helmholtz problem, both of these features are 
slowed considerably as they move into the central regions 
of the cloud, where the gas density becomes large. 

Unlike the edge-on case, the third major feature is not a 
rarefaction wave, but rather a reflected shock, which gains 
in strength as the inward propagating shock moves into 
the increasingly dense material. Thus this feature initially 
moves leftward, at a physical velocity of ~ —100 km/s, but 
a positive velocity ~ 100 km/s in the frame of the incom- 
ing shock. By f = 5 however, this shock attains a positive 
physical velocity 80 km/s, and exhibits a density ratio 
~ 4. 

In the lower three panels of these figures, we again find 

an increase in the turbulent kinetic energy along the in- 
terface, but with a lower overall magnitude than in the 
shear layer case. In fact, the most efficient source of K 
in this simulation is the backwards moving shock, which 
is made up completely of enriched gas. This means that 
mixing in this geometry is less efficient than in the edge-on 
case, and our results indicate that only ^ ^ 3% of the total 
cloud mass is able to be enriched before our final time step. 
Comparing this mixing to an identical run in which Vt was 
set to zero in eq. (39), we find that much of this mixing is 
due to numerical diffusion. At least in the one-dimensional 
case, mixing is dominated by shear, and occurs at a level 
^ ^ 20% before the onset of star formation. 

From these calculations we expect a reasonable value for 
the fraction of impinging metals mixed throughly into the 
cloud to be ^ ~ 50%. Furthermore, even higher ^ values 
may be appropriate when the motion of the ffow is not 
constricted to such limited geometries, as hinted at by a 
feature in our face-on calculations. In the regions between 
the material interface and the reverse shock, a large tur- 
bulence scale length is generated hy t ^ 1 and grows with 
time, suggesting that a global instability may be develop- 
ing within the cloud. In fact, while not strictly turbulence, 
a global mean-flow distortion is known to arise in an anal- 
ogous configuration in the laboratory, where it serves as a 
prime example of the problem of transition to turbulence. 
In the experimental studies described in Vorobieff et al. 
(2003), a column of heavy gas (SFg) initially surrounded 
by air was accelerated by a planar shock. In this case the 
cloud is observed to deform as an arc, which buckles along 
the edges. Eventually these edges spiral over themselves, 
leading to two large counter-rotating vortices as shown in 
Figure 11. These vortices quickly transition to an overall 
turbulent flow, unlike the one seen in our one- dimensional 
studies. 

It is this global distortion that is likely to be key in in- 
creasing mixing in the post-shock cloud. Indeed, a similar 
geometry was hinted at by the arcs observed in the jet- 
cloud interactions studied by Fragile et al. (2004), as well 
as by the the large-scale vortices observed in the shock- 
cloud simulations of Klein, McKee, & Colella (1994). How- 
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ever the properties of this overall buckling will undoubt- 
cidly be modified by the presence of a radial density gradi- 
ent, forces due to the dark-matter potential, and self grav- 
ity, none of which were included in these studies. Clearly 
then, for the moment, we can only conclude that the simple 
Kelvin-Hclmholtz and Rayleigh- Taylor instabilities simu- 
lated in our one-dimensional calculations are insufficient 
to produce the observed homogeneity within globular clus- 
ters. The question of global mixing requires more exten- 
sive modeling, an issue we are actively pursuing and intend 
to describe in a future publication. 

5. DISCUSSION 

The first cosmic structures to form are likely not to have 
been galaxies, but rather lower mass clouds of gas and dark 
matter too small to cool by atomic processes. In fact, in 
the Cold Dark Matter cosmological model assumed here, 
as much as 15% of the total mass is contained in such 
minihalos by a redshift of 10. Thus these clouds represent 
a vast reservoir of material that can quickly be converted 
into stars though shocking and enrichment by neighboring 
objects. 

A second generic feature of any hierarchical model is a 
high-redshift epoch of galaxy outfiows. Although the sizes 
of galaxies increase strongly with time, internal properties 
such as the scale of OB associations (McKee & Williams 
1997) and the efficiency of energy deposition from super- 
novae are largely independent of redshift. Thus while the 
ejecta of conglomerations of Type II SNe are primarily 
confined to the interstellar medium of large galaxies like 
the Milky Way forming at z ~ 2, they are easily able to es- 
cape from the dwarf galaxies forming at z ^ 3, generating 
large outfiows of material such as are now well-observed in 
the Lyman break population. 

It is clear that these outflows will interact with miniha- 
los, the uncertainty lies only in the nature of that inter- 
action. Through our estimates in §3, we have shown that 
for a large range of outflow energies, minihalos sizes, and 
separations, the most likely outcome is a dense cluster of 
stars that is bounded by self-gravity, but stripped from its 
associated dark matter. Furthermore the present-day imi- 
verse has natural candidates for these collections of stars, 
which formed on the periphery of galaxies. 

Although globular cluster formation may well be an on- 
going process, the majority of halo globular star clusters 
are old, and exhibit several features of induced minihalo 
star-formation. Observations of GC tidal streams show no 
indication of associated dark mater halos, unlike the ac- 
creted dwarf spheroidal galaxies {e.g. Ibata & Lewis 1998; 
Mayer et al. 2002). GCs also exhibit a maximum mass 
^ 10*' Mq that can not be understood by dynamical fric- 
tion or any other known destruction mechanism that oper- 
ates in the Galaxy. Instead, this value closely corresponds 
to the minimum gas mass of high-redshift minihalos for 
which atomic cooling at the virial temperature is effective. 
Finally, the span of time At over which halo GCs form in 
our model is fairly short, on the order the parent galaxy's 
virial radius divided by its outflow velocity. Even for a 
large galaxy like the Milky Way, At ^ 1 Gyr, which is com- 
parable to the observed age spread. Additionally, Freeman 
& Bland-Hawthorn (2002) have remarked on the fact that 
halo GCs and stars in the thick disk of our Galaxy have 
very similar abundances and ages. While some of these 



disk stars may be from globulars long since destroyed, ac- 
cording to our picture, others may have formed from gas 
in outflow-enriched protocluster clouds for which the con- 
ditions were never conducive to GC formation. 

The small range in metallicitics amongst the stars within 
any given cluster points to an external mechanism for en- 
richment, as provided in our picture. Furthermore the 
time scale for this enrichment must be short, as the pres- 
ence of dust and metals will inevitably spark cooling and 
star formation, which will further enrich material that is 
not quickly formed into stars. Outflow-minihalo interac- 
tions have the potential to achieve such rapid mixing, but 
order of magnitude estimates can not confirm this defini- 
tively. Rather the problem of mixing is a delicate one, in 
which the time scale for star formation in shocked regions 
is comparable to the global dynamical time, and the de- 
tailed turbulent structure is key. In the numerical studies 
presented in this paper we have taken a first step towards 
improving these estimates, exploring a simple model of 
turbulent mixing in a one-dimensional context. Yet the 
^ ^ 20% shear-induced mixing seen in this study high- 
lights the importance of large scale vortices in turbulence 
generation. Multi-dimensional simulations are necessary 
to move further, and a definitive answer will require de- 
tailed modeling of cooling, gravitation, and perhaps even 
the turbulent spectral distribution. 

Although the issue of mixing is complicated, two prop- 
erties of the metals available in outflow-minihalo interac- 
tions are suggestive of recent GC observations. First, as 
pointed out by Beasley et al. (2003), no single Pop III 
object seems capable of producing the typical abundance 
ratios of heavy elements in GC stars. Some distribution 
of flrst generation stars is needed, and a galactic outflow 
is the likely consequence of several concurrent supernovae. 
Second, the existence of a correlation between a parent 
galaxy's luminosity and the mean color (viz., metallicity) 
of its GC systems (Strader, Brodie, & Forbes 2004) ar- 
gues strongly against self-enrichment as the source of a 
GCs heavy elements. 

Our basic model even enables us to deduce such a de- 
pendence. As the efficiency of our scenario is primarily 
due to an interplay between cooling and sound-crossing 
times, both which depend on the shock velocity, we fix 
a typical at which GC formation is effective. From 
eq. (1) this gives that Rs cx '^(1 + Zg)^^. We then as- 
sume that, on average, the total mass in metals and the 
outflow energy will be proportional to luminosity of the 
parent galaxy to which the resulting globular cluster be- 
longs, yielding Zc oc L^^f (l-l-Zs)^. From the simple top-hat 

collapse model (1 -|- Zs) oc (Mgai)~^"+^)/^, where n is the 
slope of the primordial power spectrum, so again taking 
the mass proportional to igai, we obtain the prediction 

Zc (X L-Jr"^/\ (44) 

where for the CDM model n is between -3 and and -2 
for the scales of interest. This gives a final index ranging 

from 1/3 for small galaxies to for the largest ones, brack- 
eting the observed value of 1/6 (Strader, Brodie, & Forbes 
2004). 

A similar calculation can be carried out to estimate the 
maximum galactocentric radius of the distribution of halo 
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GCs. Fixing wc have 

i?.,max oc E'/^l + zs)-' oc L^^+^^Z^ (45) 

such that for the n values appropriate for CDM, the index 
ranges from 1/3 for small galaxies to 1/2 for large galax- 
ies, and the "best-fit" n value to the observed metallicity- 
luminosity relationship gives i?s,max oc L^/^"^. As such ob- 
servations require both high-angular resolution and a wide 
field of view, no current constraints on this relationship are 
available in the literature. However, such measurements 
are well within the capabilities of the Advanced Camera for 
Surveys on the Hubble Space Telescope and new informa- 
tion in this area should be expected to become available 
soon. Note that this scaling with luminosity is likely to 
persist even if substantial "shuffling" of globular clusters 
within this maximum radius is likely to have taken place, 
as a large fraction of halo GCs are on close to radial orbits 
(eg. Dinescu et al. 1999; Wang & Zhou 2003). On the other 
hand, any initial metallicity gradient established among a 
galaxy's halo GCs is expected to be washed-out by this 
shuffling and possibly by the addition of clusters striped 
from small, neighboring galaxies. It is also interesting to 
note that while no clear metallicity gradients have ever 
been observed in halo GCs, there are hints of a fall off in 
metallicity with radius in the Milky Way Halo GC popula- 
tion when one attempts to remove accreted systems (Zinn 
1993; Parmentier et al. 2001). 
While these comparisons are perhaps a bit simplified. 



it is reassuring that recent investigations not only uncov- 
ered a trend with luminosity, but one compatible with our 
model, and that we are able to construct a simple pre- 
diction that will be testable in the near future. Perhaps 
even more reassuring is the fact that the ideas presented 
here allow for immediate points of contact between theory 
and new observations. Although many issues remain to be 
better understood, our scenario provides clear directions 
that can lead to more detailed comparisons with various 
observational constraints. Such future studies may soon 
shed new light on the old question of the origin of halo 
globular clusters. 

We are grateful to Jean Brodie, Alex Heger, Chris 
McKee, Francesco Palla, Jay Stradcr, and an anony- 
mous referee for helpful comments during the prepara- 
tion of this manuscript. Thanks with much appreciation 
to Christopher D. Tomkins, Hydrodynamics Group DX-3, 
Los Alamos National Laboratory, who provided us with 
the experimental results presented in Figure 11. ES was 
supported by an NSF Math and Physical Sciences Distin- 
guished International Postdoctoral Research (NFS MPS- 
DRF) fellowship during part of this investigation; his re- 
search was also supported by the National Science Foun- 
dation under grant PHY99-07949. The work performed by 
JW and FH was conducted under the auspices of the U. 
S. Department of Energy's contract W-7405-ENG-36 with 
the University of California. 



REFERENCES 



Abel, T., Bryan, G. L., & Norman, M. L. 2002, Sci, 295, 93 
Addessio, F. L., Baumgardner, J. R., Dukowicz, J. K., Johnson, 
N. L., Kashiwa, B. A., Raucnzahn, R. M., & Zcmach, C. 1990, 
CAVEAT: A Computer Code for Fluid Dynamics Problems with 
Large Distortion and Internal Slip, LA-MS-10613, Los Alamos 
National Laboratory 
Armandroff, T. C, 1989, AJ, 97 375 
Ashman, K. M., & Bird, C. M. 1993, AJ, 106, 2281 
Ashman, K. M., & Zepf, S. E. 1998, Globular Cluster Systems, (New 

York : Cambridge University Press) 
Beasley, M. A., et al. 2000, MNRAS, 318, 1249 

Beasley, M. A., Kawata, D., Pcarce, F. R., Forbes, D. A., & Gibson, 

B. K. 2003, ApJ, 596, L187 
Besnard, D., Harlow, F. H., Rauenzahn, R. M., & Zemach, C. 1992, 

Turbulence Transport Equations for Variable-Density Turbulence 

and Their Relationship to Two-Field Models, LA-MS-12303, Los 

Alamos National Laboratory 
Blais-Ouellette, S., Carignan, C, Amram, P., & Cote, S. 1999, AJ, 

118, 2123 

Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., Kravtsov, 
A. v.. Klypin, A. A., Primack, J. R., & Dekel, A. 2001, MNRAS, 
321, 559 

Bromm, V., Coppi, P. S., Larson, R. B. 2002, ApJ, 564, 23 

Bromm, V., & Clarke, C. J. 2002, ApJ, 556, LI 

Brown, J. H., Burkert, A., & Truran, J. W. 1995, ApJ, 440, 666 

Caputo, F., & Castellani, V. 1984, MNRAS,207, 185 

Carignan, C, & Beaulieu, S. 1989, ApJ, 347, 760 

Carignan, C, Sancisi, R., & van Albada, T. S. 1988, AJ, 95, 37 

Cen, R. 2001, ApJ, 560, 592 

Ciardi, B., Ferrara, A., Governato, F., & Jenkins, A. 2000, MNRAS, 
314, 611 

Dalganro, A., & McCray, R. A. 1972, ARA&A, 10, 375 
de Blok, W. J. G., & McGaugh, S. S. 1997, MNRAS, 290, 533 
Dinescu, D. I., Girard, T. M., van Altena,W. F. 1999, ApJ, 117, 
1792 

Eke, V. R., Cole, S., & Frenk, C. S. 1996, MNRAS, 282, 263 

Eisenstein, D. & Hu, W. 1999, ApJ, 511, 5 

Elmegreen, B. G., & Efremov, Y. N. 1997, ApJ. 480, 235 

Fall. M. S., & Rees, M. J. 1985, ApJ, 298, 18 

Ferrara, A. 1998, ApJ, 499, LI 

Forbes,'D. A., Brodie,' J. P., & Huchra, J. P. 1997, AJ, 113, 887 



Fragile, P. C, Murray, S. D.. Anninos, P., & van Breugel, W. 2004, 
ApJ, 604, 74 

Freeman, K. & Bland-Hawthorn, J. 2002, ARA&A, 40, 487 
Fryc, B., Broadhurst, T., & Benitez, N. 2002, ApJ, 568, 558 
Fujita, A., Martin, C. L., Mac Low, M.-M., & Abel. T. 2003, ApJ, 
599. 50 

Fujita, A., Mac Low, M.-M., Ferrara, A., & Meiksin, A. 2004, ApJ, 

accepted (astro-ph/0405611) 
Galli, D. & Palla, F. 1998, A&A, 335, 403 
Gnedin, N. Y. 2000 ApJ, 535 530 

Grillmair, C. J., Freeman, K. C, Irwin, M. & Quinn, P. J. 1995, AJ, 

109, 2553 

Gunn, J. E. 1980 in Globular Clusters, eds. D.Hanes and B. Madore 

(Cambridge University Press), 301 
Haiman, Z., Rees, M., & Loeb, A. 1997, ApJ, 476, 458 (erratum 

484, 985) 

Harlow, F. H., & Amsden, A. A. 1971, Fluid Dynamics a LASL 

Monograph LA-MS-4700, Los Alamos Scientific Laboratory 
Harlow, F.H., & Nakayama, P.I. 1967, Phys. Fluids, 10, 2323 
Heger, A. & Woosley, S. E. 2002, ApJ, 567, 532 
Hu, E. M., et al. 2002, ApJ, 568, L75 

Ibata, R., Irwin, M., Lewis, G. F., & Stolte, A. 2001, ApJ, 547 L133 

Ibata, R. A. & Lewis, G. F. 1998, ApJ, 500, 575 

Irwin, M. J.& Hatzidimitriou, D. 1993 in ASP Conf. Ser. 48, The 

Globular Cluster-Galaxy Connection, eds. G. H. Smith & J. P. 

Brodie (San Francisco: ASP), 322 
Kitayama, T., Yoshida, N., Susa, H., & Umemura, M. 2004 ApJ, in 

press, (astro-ph/0405280) 
Klein, R. , McKce, C, & Colella, P. 1994, 420, 213 
Kogut, A. et al. 2003, ApJS, 148, 161 
Krauss, L. M., & Chaboycr, B. 2003, Science, 299, 65 
Lacey, C. & Cole, S., 1993, MNRAS, 262, 627 

Larsen, S. S., Brodie, J. P., Huchra, J. P., Forbes, D. A., Grillmair, 

C. J. 2001, AJ, 121, 2974 
Launder, B. E., Reece, G. J., & Rodi, W. 1975, J. Fluid Mech, 68, 

537 

Li, Y., & Burstein, D., 2003, ApJ, 598, L103 
Mac Low, M.-M., & Slmll, J. M. 1986, ApJ, 302, 585 
Mac Low, M.-M., & Ferrara, A. 1999, ApJ, 513, 142 
Madau, P., Ferrara, A., Rees, M. J. 2001, ApJ, 555, 92 
Marm-Franch, A., & Aparicio, A. 2003, ApJ, 585, 714 
Marinoni, C, & Hudson, M. J. 2002, ApJ, 569, 101 



13 



Mashchenko, S. & Sills, A. ApJ, submitted 

Mayer, L., Moore, B., Quinn, T., Governato, F., &; Stadel, J. 2002, 

MNRAS, 336, 119 
McKee, C. F., & Williams, J. P. 1997, ApJ, 476, 144 
McLaughlin, D. E. 1999, AJ, 117, 2398 
Moore, B. 1996, ApJ, 461, LIS 

Mori, M., Fcrrara, A., & Madau P. 2002, ApJ, 571, 40 
Murray, S. D., & Lin, D. N. C. 1992, ApJ, 400, 265 
Nakasato, N., Mori, M., & Nomoto, K. 2000, ApJ, 535, 776 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 
(NFW) 

Ostriker, J. P. , Spitzer, L. & Chevalier, R. A. 1972, ApJ, 176, L51 

Ott, J., Martin, C. L., & Walter, F. 2003, ApJ, 594, 776 

Palla, F. & Zinnecker, H. 1988, in lAU Symposium 126: The Harlow 

Shapley Symposium on Globular Cluster Systems in Galaxies, eds. 

J. E. Grindlay & A.G. D, Philip (DoredrechtiRcidcl), p323 
Parmenticr, G., Jehin, E., Magain, P., Noels, A., & Thoul, A. A. 

2000, A &A, 363, 526 
Peebles, P. J. E., & Dicke, R. H. 1968, ApJ, 154, 891 
Peng, W. & Weisheit, J. C. 1991, PASP, 103, 891 
Perlmutter, S. et al. 1999, ApJ, 517, 565 
Peterson, J. R. ct al. 2003, ApJ, 590, 207 

Pettini, M., Stcidcl, C. C, Adclbcrger, K. L., Dickinson, M., & 
Giavalisco, M. 2001, ApJ, 528, 96 

Press, W. H., & Schcchtcr, P. 1974, ApJ, 187, 425 (PS) 

Rees, M. 2003, in Star and Structure Formation: From First 
Light to the Milky Way. (available at http://www.exp- 
astro.phys . cthz . ch/ETH_Astro_2003/Talks /reesj .pdf ) 

Renzini, A. 1997, ApJ, 488, 35 

Scannapicco, E., Fcrrara, A., & Broadhurst, T. 2000, ApJ, 536, Lll 
Scannapicco, E., Fcrrara, A., & Madau, P. 2002, ApJ, 574, 590 
Scannapicco, E., Schneider, R., & Fcrrara, A. 2003, ApJ, 589, 35 
Schwcizcr, F., Miller, B. W., Whitmorc, B. C, & Fall, S. M. 1996, 

AJ,112, 1839 
Searle, L., & Zinn, R. 1978, ApJ, 225, 357 
Shapiro, P., R., & Kang, H. 1987, ApJ, 318, 32 
Shapiro, P., R., Iliev, L T., & Raga, A. C. 2004, MNRAS, 348, 753 
Songaila, A. 2001, ApJ, 561, L153 (801) 



Spergel, D. N. et al. 2003, ApJS, 148, 175 
Spitzer, L., & Thuan, T. X. 1972, ApJ, 175, 31 
Strader, J., Brodie, J. P., & Forbes, D. A. 2004, AJ, in press (astro- 
ph/0403160) 

Steinkamp, M., Clark, T., & Harlow F. (1999a), Int. J. of Multiphase 

Flow, 25, 599 

Steinkamp, M., Clark, T., & Harlow F. (1999b), Int. J. of Multiphase 
Flow, 25, 639 

Suntzeff, N. 1993, in ASP Conf. Ser. 48, The Globular Cluster- 
Galaxy Connection, ed. G. H. Smith & J. P. Brodie (San Francisco: 
ASP), 167 

Sutherland, R. S., & Dopita, M. A. 1993, ApJS, 88, 253 
Swaters, R. A., Madore, B. F., & Trewhella, M. 2000, ApJ, 531, 
L107 

Swaters, R. A., Madore, B. F., van den Bosch, F. C, & Balcells, M. 

2003, ApJ, 583, 732 
Tegmark, M., Silk, J., & Evrard, A. 1993, ApJ, 417, 54 
Thacker, R. J., Scannapicco, E., & Davis, M. 2002, ApJ, 581, 836 
Tytler, D., Fan, X. M., Buries, S. et al., 1995, Q30 Absorption Lines, 
Proceedings of the ESQ Workshop Held at Garching, Germany, 21 
- 24 November 1994, edited by Georges Meylan. Springer- Verlag 
Berlin Heidelberg New York. Also ESQ Astrophysics Symposia, p. 
289 

Ranch, M., Sargent, W. L., Womblc, D. S., & Barlow, T. A. 1996, 

ApJ, 467, L5 
Uehara, H., & Inutsuka, S. 2000, ApJ, 531, L91 
van Breugel, W., Filippenko, A. V., Heckman, T., & Miley, G. 1985, 

ApJ, 293, 83 

van den Bosch, F. C, & Swaters, R. A. 2001, MNRAS, 325, 1017 
Vorobieff, P., Mohamed, N.-G., Tomkins, C, Goddcnough, C, 
Marr-Lyon, M., & Benjamin, R. F. 2003, Phys. Rev. E 68, 65301 
Wang, L. & Zhou, H.-N. 2003, CliA&A, 27, 42 
Whitmorc, B. C, & Schwcizcr, F. 1995, AJ, 109, 96 
Wooslcy S. E. & Weaver, T. A. 1995, ApJS, 101, 81 
Yoon, S.-J., & Lee, Y.-W. 2002, 297, 578 
Zinn, R. 1985, ApJ, 293, 424 

Zinn, R. 1993, Stellar Populations, eds. C. A. Norman, A. Renzini, 
& M. Tosi (Cambridge: Cambridge University Press, 73 



14 




0.2 0.4 0.6 0.8 1 

log(b/Rj 

Fig. 1.— Fiducial model of globular cluster formation (E55 = 10, (5 = 44, e = 0.3, M = log{Zc/ZQ) = -1.5). Top panel: 

fractional gas mass (M(> 6) = M(> 6)/Mgas) outside a given impact parameter. Second panel: Local metallicity as function of b. Third 
panel: Cloud momentum surface density, pc{b), (solid line) versus impinging momentum per unit area (dotted line), both plotted in units of 
Mq kms~^ kpc"'^. Fourth panel: Comparison of the sound-crossing time (dotted), free-fall time (solid), and cooling time (dashed) just after 
the shock has passed over the minihalo. Fifth panel: Comparison of the sound-crossing time, free-fall time, and cooling time after the cloud 
has cooled to lOOK. Lines are as in the fourth panel, while the long-dashed line is the cooling time at lOOOK. 
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Fig. 2. — Effect of varying the distance from the outflowing galaxy. Distances are chosen such that the overafl mass average metallicities 
are log{Zc/ZQ) = -3.0, -2.5, -2, -1.5, and -1. Curves are as in the top four panels of Fig. 1. As metallicity increases, the shock reaches 
the minihalos with more momentum and energy, increasing both ps (dotted lines, center panel) and tcool (dashed lines, bottom panel). As 
the cloud is shocked to higher temnperatures, increasing Zc decreases the sound crossing time. Minihalo clouds are disrupted in the highest 
metallicity case, while shocking is too weak to trigger star formation in the lowest-metallicity case. 
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Fig. 3. — Effect of varying the minihalo mass. Panels are as in Fig. 2, in all cases curves are for masses of and 
Increasing the mass while fixing Zc raises pc (solid lines, center panel) dramatically and moves the minihalo closer, increasing the impinging 
momentum surface density to a lesser extent (dotted lines, center panel). Similarly increasing the minihalo mass increases the cooling time 
dramatically (dashed lines, lower panel) and although the cloud is hotter, the increase in mass causes the sound crossing time to increase 
slightly. Thus shock induced star formation is most efficient in lower-mass clouds, while long cooling times suppress star formation in the 
largest minihalos. 
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Fig. 4. — Effect of varying the energy input {-B55) and tfie wind efficiency (e). Panels as in Fig. 3. Left Column: Er^^ is varied between 5 and 
30, keeping Zc fixed. The largest value of £55 corresponds to the largest impinging momentum surface density in the center panel (shaded 
region), but to the lowest post-shock temperature. Thus E55 = 30 corresponds to the largest value of tsc (shaded region bottom panel) and 
the smallest value of tcool (cross-hatched region bottom panel). This is because objects of a fixed metallicity lie at larger distances as £55 
is increased. Right Column: e is varied between .1 and .5. As with E^s, an increase in the efficiency parameter increases ps (shaded region, 
center panel), but in this case tsc decreases (shaded region bottom panel), and tcool goes up (cross-hatched region bottom panel). 
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Fig. 5. — Effect of varying the density surrounding the outflowing galaxy (S) and the cloud virialization redshift (zc), panel and curves as 
in Fig. 4. Left Column: S is varied between 10 and 80. Increasing S increases ps and reduces the post-shock temperature. Thus as 5 goes up, 
*cool decreases (cross hatched region, lower panel), and tuc goes up (shaded region, lower panel), improving the formation efficiency. Right 
Column: Zc is varied between 8 and 15. This has essentially the opposite effect from 5. Thus increasing Zc decreases ps, increases tcooli ^.nd 
reduces tac, making our mechanism less efficient. Finally, increaing Zc also results in a more compact cloud, with a higher escape velocity 
(solid region, center panel) and a shorter free-fall time (solid region, bottom panel). 




Fig. 6. — Illustration of the edge-on (left panel) and face-on (right panel) configurations. In both panels the cloud of primordial composition 
is indicated in black, while the impinging enriched material is indicated in grey. 
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Fig. 7. — Quantities in the Kelvin-Helmholtz shear problem at times t = (left column) and t = 1 (1.5 X 10^ years, right column). In the 
top row the solid line is c\, the long-dashed line is C2, and the short-dashed line is a measure of the normalized cloud mass M{b) exterior to 
an impact parameter, b. In the second row, the units of p are grams per cm^. In the third through fifth rows Ux, Uy, and Cs are in km/s. 
The last three rows represent turbulent quantities. Among these K (sixth row), is in ergs/gram, S (seventh row) is in units of parsecs, and 
ax and ay (represented by the solid and dashed lines in the final row) are in km/s. The initial turbulence kinetic energy is assumed to bo 2% 
of initial total energy. 
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Fig. 8. — Kelvin-Helmholtz shear problem at times t = 2 (3.0 X 10® yr, left column) and t = 5 (7.5 X lO'' years, right column). Rows are 
as in Figure 7, and the additional dotted lines in the uppermost t = 5 panel show ci and C2 in a case in which ut has been set to zero in the 
transport equations for the concentrations, and the lower dashed curve shows M[i) for this case. 
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Fig. 9. — Face-on Raleigh- Taylor problem at times t = (left column) and t = \ (1.5 X 10® years, right column). In the upper row the solid 
line is Ci, the long-dashed line is C2, and the short-dashed line is the normalized cloud mass MiC) contained in front of a given distance from 
the back of the cloud, t. Rows 2-4 represent mean flow quantities, and rows 5-7 represent turbulent quantities, with units as is 7. 
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Fig. 10. — Face-on Raleigh-Taylor problem at times t = 1 (3.0 X 10® yr, left column) and t = 5 (7.5 X 10® yr, right column). Rows are as 
in Figure 9, and, as in Figure 8, additional dotted lines have been added in the upper right panel which show C\ and C2 in a case in which 
vt has been set to zero in the transport equations for the concentrations. 
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Fig. 11. — Experimental results of a column of heavy gas (SFg) which is overtaken by a planar shock, which moves in from the left in 
these figures. The column is 3.1 mm in diameter and time frames axe at 320 (left) and 470 (right) microseconds. In our units of radius over 
the shock speed, these correspond to times of 90 and 130 repectively. This image was provided by C. D. Tomkins and was taken from the 
experiments described in Vorobieff et al. (2003). 



